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Along the line of thoughts of Berry and Robnik,!^ we investigated the gap distribution 
function of systems with infinitely many independent components, and discussed the level- 
spacing distribution of classically integrable quantum systems. The level spacing distribution 
is classified into three cases: Case 1: Poissonian if ju(+oo) = 0, Case 2: Poissonian for large 
S, but possibly not for small S if < p,(+oo) < 1, and Case 3: sub-Poissonian if p,(+oo) = 1. 
Thus, even when the energy levels of individual components are statistically independent, 
non-Poisson level spacing distributions are possible. 

An important property of quantum-classical correspondence appears in the sta- 
tistical property of energy levels of bounded quantum systems in the semiclassical 
limit. Universal behaviors are found in the statistics of unfolded energy levels at 
a given interval, which are the sequence of numbers uniquely determined by the 
energy levels using the mean level density obtained from the Thomas-Fermi rule. 
It is widely known that, for quantum systems whose classical counterparts are in- 
tegrable, the distribution of nearest-neighbor level spacing is characterized by the 
Poisson distribution!), while for quantum systems whose classical counterparts are 
strongly chaotic, the level statistics are well characterized by the random matrix 
theory which gives level-spacing distribution obeying the Wigner distribution. eP 

Level statistics for the integrable quantum systems has been theoretically stud- 
ied by Berry-TaborcP, SinaiuP, MolchanovEP, BlehertP, Connors and Keating©, and 
MarkloftP, and have been the subject of many numerical investigations. Still, its 
mechanism is not well understood, the appearance of the Poisson distributions is 
now widely admitted as a universal phenomenon in generic integrable quantum sys- 
tems. 

As suggested, e.g., by Hannay (see the discussion of Ref.l)), one possible expla- 
nation would be as follows: For an integrable system of / degrees-of-freedom, almost 
every orbit is generically confined in each inherent torus, and the whole region in 
the phase space is densely covered by invariant tori as suggested by the Liouville- 
Arnold theorem. e$ In other words, the phase space of the integrable system consists 
of infinitely many tori which have infinitesimal volumes in Liouville measure. The 
energy level sequence of the whole system is then a superposition of sub-sequences 
which are contributed from those regions. Therefore, if the mean level spacing of 
each independent subset is large, one would expect the Poisson distribution as a 
result of the law of small numbers. c3 This scenario suggested by Hannay is based 
on the theory proposed by Berry and Robnik.&P 



typeset using VTpTfX.sty <ver.l.0> 



2 



H. Makino and S. Tasaki 



The Berry-Robnik theory relates the statistics of the energy level distribution 
to the phase-space geometry by assuming that the sequence of the energy spectrum 
is given by the superposition of statistically independent subspectra, which are con- 
tributed respectively from eigenfunctions localized onto the invariant regions in phase 
space. Formation of such independent subspectra is a consequence of the condensa- 
tion of energy eigenfunctions on disjoint regions in the classical phase space and of 
the lack of mutual overlap between their eigenfunctions, and, thus, can be expected 
only in the semi-classical limit where the Planck constant tends to zero, Ti — > 0. 
This mechanism is sometimes referred to as the principle of uniform semi- classical 
condensation of eigenstates^^ which is based on an implicit state by Berry. 

In this paper, keeping the above mentioned scenario in mind, we derive the gap 
distribution function of systems with infinitely many components, and discuss the 
level spacing statistics of integrable quantum systems. 

In the Berry-Robnik approach, IIP the overall level spacing distribution is derived 
as follows: Consider a system whose classical phase space is decomposed into N- 
disjoint regions. The Liouville measures of these regions are denoted by pi(i = 
1, 2, 3, ■ ■ ■ , N) which satisfy J2iLi Pi = 1- Let E(S) be the gap distribution which 
stands for the probability that an interval (0, S) contains no level. E(S) is expressed 
by the level spacing distribution P(S) as, E(S) = da f£° P(x)dx. When the 
entire sequence of energy levels is a product of statistically independent superposition 
of N sub-sequences, E(S; N) is decomposed into those of sub-sequences, Ei(S; pi), 

N 

E(S;N) = l[E i (S;p i ). (0-1) 

i=l 

In terms of the normalized level spacing distribution pi(S;pi) of a sub-sequence, 
Ei(S;pj) is given by Ei(S;pi) = pi da f™ pi(x; Pi)dx, and pi(S;pi) is assumed to 
satisfy EP 

S- Pi (S; Pi )dS = -. (0-2) 
Jo Pi 

This equation is satisfactory when the Thomas- Fermi rule for individual phase space 
regions still holds. 

Note that the spectral components are not always unfolded automatically in 
general even when the total spectrum is unfolded. However, in the sufficient small 
energy interval, each spectral component obeys a same scaling law (see Appendix A 
of Ref.14)), and thus is unfolded automatically by an overall unfolding procedure. 
Eqs. (|0-1| ) and (|0-2|) relate the level statistics in the semiclassical limit with the 
phase-space geometry. 

In most general cases, the level spacing distribution might be singular. In 
such a case, it is convenient to use its cumulative distribution functions: Pi(S) = 
Jo Pi{x;pi)dx. 

In addition to Eqs ( pT| ) and Q0-2| ), we assume the following two conditions: 

• Assumption (i): The statistical weights of independent regions uniformly van- 
ishes in the limit of infinitely many regions: maxj pi — > as ./V — > +oo. 

• Assumption (ii): The weighted mean of the cumulative distribution of energy 
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+oo to p(p). The limit is 



spacing, p(p;N) = T,i=i PiP>i(p) , converges in N 

uniform on each closed interval: < p < S. 
Under assumptions (i) and (ii), eqs.(OT) and ( |0-2j ) lead to the overall level spacing 
distribution whose gap distribution function is given by the following formula in the 
limit of N — > +00, 

r s 

E fl (S)=exp - (1 -fl(a)) da , (0-3) 
Jo 

where the convergence is in the sense of the weak limit. When the level spacing dis- 
tributions of individual components are sparse enough, one may expect p = and 
the gap distribution of the whole energy sequence is reduced to the Poisson distri- 
bution, Ep=o(S) = exp (-S). In general, one may expect p 7^ which corresponds 
to a certain accumulation of the levels of individual components. 

In what follows, starting from eqs.([H|) and (p-2|), and assumptions (i) and (ii), 
Eq.( |0-3| ) is derived in the limit of N — ► +00, and by analyzing Eq.( p-3[ ), the level 
spacing distribution is discussed. 

Following the procedure by Mehta,© we rewrite E(S; N) in terms of the cumu- 
lative level spacing distribution pi(S) of independent components: 



N r 

E(S;N) = H pi / 

i=l L Js 



Pi l da (I - pi(a)) 
s 



N 

n 

i=l 



Pi 



da(l - Pi(a)) 



(04) 



The second equality follows from Eq. (p-2|) , integration by parts and hnv^+oo a (1 — Pi(a)) 
0, which results from the existence of the average. Since the convergence of J2iLi Pi^i{v) — ¥ 
p(a) for N — > +00 is uniform on each interval a E [0,5] by Assumption (ii), and 
|^i(o")| < 1) E(S; N) has the following limit in N — > +00: 



s ^ s 

log E(S;N) = - I da[l-p,(a;N)]+J20(p?) — . - / da [1 
Jo i Jo 



e + 0(e 2 ) in e <C 1, and the follow- 



where we applied the expansion log(l + e) = 

ing property obtained from Assumption (i): | J2iL\ 0(pf)\ < C ■ maxj p,i J2iLi Pi 
C ■ maxj pi — > as N — > +00 with C a positive constant. Therefore, we have 
Eq.(|0-3|). We remark that, when p(S) is differentiable, the asymptotic level spacing 



distribution is described as Pp,(S) = [(1 — p(S)) 2 + p'(S)] exp — J Q (1 — p(a)) da 

Since Pi(S) is monotonically increasing and < Pi(S) < 1, p(S) has the same 
properties. Then, 1 — p(S) > for any S > and one has 4 Jo da (I — p{a)) — > 
1 — ^(+00) as S — > +00. According to this limit, the level spacing distribution 
corresponding to eq. (|0-3| ) is classified into the following three cases: 



Case 1, /2(+oo) = 0: The limiting level spacing distribution is the Poisson 
distribution. Note that this condition is equivalent to p(S) = for v 5 because 
p(S) is monotonically increasing. 

Case 2, < p(+oo) < 1: For large S values, the limiting level spacing distribu- 
tion is well approximated by the Poisson distribution, while, for small S values, 
it may deviate from the Poisson distribution. 
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• Case 3, /2(+oo) = 1: The limiting level spacing distribution deviates from the 
Poisson distribution for V S', and decays as S — > +oo more slowly than does the 
Poisson distribution. This case will be referred to as a sub-Poisson distribution. 
One has Case 1 if the individual level spacing distributions are derived from scaled 
distribution functions fa as pi(S;pi) = pifi(piS), where /, satisfy J +o ° fi(x)dx = 
1 and / +O ° xfi(x)dx = 1, and are uniformly bounded by a positive constant D: 
\fi{S)\ < D (1 < i < N and S > 0). Indeed, one then has 

TV s N 

HS;N)\<y2p 2 t / |/ i ( ft x)| ( ix< J D5maxp i Vp J ^0 = /i(5). 

In general, one may expect ji(S) ^ which corresponds to the non-Poisson distribu- 
tion. Such a case is expected when there is strong accumulation of the energy levels 
of individual components which leads to the non-smooth cumulative distribution 
function pi(S). For a certain system class^such accumulation is observable. One 
known example is the rectangular billiard. EJ) The level spacing distribution of this 
system shows strong accumulation when the aspect ratio of two sides of a billiard wall 
is close to a rational. Another example is the two-dimensional harmonic oscillator 
whose level spacing distribution is non-smooth for arbitrary system parameter, a The 
final example is studied by Shnirelman, Chirikov and Shepelyansky, and Frahm and 
Shepelyansky for a certain type of system wjiich contains a quasi-degeneracy result 
from inherent symmetry(time reversibility). EZrlIs''[Lj) As is well known, the existence 
of quasi-degeneracy leads to the sharp Shnirelman peak at small level spacings. Such 
phenomena will be discussed in forthcoming papers. 

In most general cases, the integral in equation (|0-3| ) converges in S <C +oo 
and then lims->+oo Ep(S) ^ 0, the limiting gap distribution Ep_(S) does not work 
accurately. In such case, however, its differentiation still work accurately in S — * +co 
limit ( see Ref.14) ), and thus the above classification (Case 1-3) holds in general. 

In this paper, we investigated the gap distribution function of systems with 
infinitely many independent components, and discussed the level-spacing statistics 
of classicaly integrable quantum systems. In the semiclassical limit, reflecting in- 
finitely fine classical phase space structures, individual eigenfunctions are expected 
to be well localized in the phase space and contribution independently to the level 
statistics. Keeping this expectation in mind, we considered a situation in which the 
system consists of infinitely many components and each of them gives an infinitesi- 
mal contribution. And by applying the arguments of Mehta, and Berry and Robnik, 
the limiting gap distribution was obtained which is described by a single monotoni- 
cally increasing function p,(S) of the level spacing S. Three cases are distinguished: 
Case 1: Poissonian if p(+oo) = 0, Case 2: Poissonian for large S, but possibly not 
for small S if < p(+oo) < 1, and Case 3: sub-Poissonian if /2(+oo) = 1. Thus, 
even when the energy levels of individual components are statistically independent, 

non-Poisson level spacing distributions are possible. 

The authors would like to thank Professor M. Robnik and A. Shudo for their 
helpful advice. The authors also thank the Yukawa Institute for Theoretical Physics 
at Kyoto University. Discussions during the YITP workshop YITP-W-02-13 on 
"Quantum chaos: Present status of theory and experiment" were useful to com- 
plete this work. This work is partly supported by a Grant-in-Aid for Scientific 



5 



Research (C) from the Japan Society for the Promotion of Science, and by the Min- 
istry of Education, Science, Sports and Culture, Grant-in- Aid for Young Scientists 
(B)(15740244,2003-2005). 



References 

[1] M.V. Berry and M. Robnik, J. Phys. A17 (1984), 2413. 

[2] M.V. Berry and M. Tabor, Proc. R. Soc. Lond. A356(1977),375. 

[3] O. Bohigas, M.J. Giannoni, and C. Schmit, Phys. Rev. Lett. 52(1984), 1. 

[4] Ya. G. Sinai, Physica A163(1990),197. 

[5] S. A. Molchanov, Commun. Math. Phys. 78(1981),429 ; S. A. Molchanov, Math. USSR 
Izvestija 12(1978),69. 

[6] P.M. Bleher, J.Stat. Phys. 61(1990), 869 ; P.M. Bleher, J.Stat. Phys. 63(1991), 261. 
[7] R. D. Connors and J. P. Keating, J. Phys. A30(1997), 1817. 
[8] J. Marklof, Comm. Math. Phys. 199(1998), 169. 

[9] V.I. Arnold and A. Avez, Ergodic Problems of Classilcal Mechaniocs (Benjamin, N.Y.,1978) 

; V.I. Arnold, Russ. Math. Surveys 18:6(1963), 85. 
[10] W. Feller, An introduction to probability theory and its applications (2nd ed., John Wiley 

& Sons, Inc., New York, 1957) 
[11] B. Li and M. Robnik, J. Phys. A27(1994),5509. 

[12] M. Robnik, Nonlinear Phenomena in Complex Systems l,nl(1998),l. 

[13] M. V. Berry,J. Phys. A10(1977),2083. 

[14] H. Makino and S. Tasaki, Phys. Rev. E67(2003), 066205. 

[15] M.L. Mehta,i?andom Matrices(2nd ed., San Diego, CA: Academic Press, 1991) App. A. 2. 

[16] M. Robnik and G. Veble, J. Phys. A31(1998),4669. 

[17] A.I. Shnirelman, Usp.Mat.Nauk. a) 30 N4 (1975), 265. 

[18] B.V. Chirikov and D.L. Shepelyansky, Phys.Rev.Lett. 74(1995), 518. 

[19] K.M. Frahm and D.L. Shepelyansky, Phys.Rev.Lett. 78(1997), 1440. 



